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1 Introduction 



The study of back reaction of gravitational waves on black hole spacetimes 
has been hampered by technical and conceptual difficulties, from the lack of 
closed-form expressions for the perturbations to subtle issues of gauge. The 
classical problem of the back reaction of a metric perturbation propagating 
out of (or into) a black hole is very interesting, but already complicated 
enough, and any progress in the classic realm would facilitate tremendously 
the analysis of quantum effects. 

One of the main difficulties in considering the back-reaction of linearized 
perturbations is the problem of gauge [1-3]. We address this issue by consid- 
ering a consistent perturbative expansion of the Einstein equations in first- 
and second-order perturbations of the black hole metric. There are gauge 
transformations related to the parametrization of the dynamical degrees of 
freedom in each order in perturbation theory, and we use these symmetries 
judiciously to cast the problem in as simple a manner as possible. The task 
of identifying physical quantities in the effective energy-momentum tensor 
(henceforth EEMT) of gravitational perturbations becomes much simpler 
that way. 

The main point of this paper is the possibility of identifying physical 
quantities in Regge- Wheeler gauge. This issue has also been raised by Gleiser 
[|] for the case of the ^ = 2, even-parity gravitational perturbation. We 
explore the connections between the flux of energy out of the black hole, 
contained in the G*^ equations, and the inner product necessary both for the 
proper normalization of the gravitational perturbations and for any attempts 
at quantizing the gravitational fluctuations. 

The outline of the paper is as follow: in Section 2 we review the formal- 
ism connected with the perturbative expansion of the Einstein equations to 
second order and deflne the EEMT of gravitational perturbations, T^^. In 
Section 3 we write the linearized perturbations and consider a scalar product 
that can flx their normalizations. In Section 4 we give details about our flrst 
order gauge choice: the Regge- Wheeler (RW) one. In Section 5 we show 
how the flux of gravitational waves computed in the RW gauge satisfles all 
the physical criteria of an energy flux both at inflnity and on the horizon. 
In Section 6 we present a perturbative expansion of the physical degrees of 
freedom (the Regge- Wheeler and Zerilli functions) that can be used in the 
expressions for the energy, then we show that it is possible to normalize these 
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gravitational waves in the Schwarzschild background. Finally, we compute 
what that energy is in the case of the odd perturbations. We conclude in 
Section 7. 

2 The General Plan 

We study back reaction by expanding the Einstein field equations to second 
order in the initial values of the metric perturbations. The metric of the 
perturbed black hole is given by the series 

9,. = + + e'g'-^J + 0{e^) , (1) 

where e is the perturbative parameter, which can be thought of as the am- 
plitude of the metric perturbation on the initial value surface. It is useful for 
the moment to regard the first order metric perturbations as gravitational 
waves propagating in the black hole background, and second order metric 
perturbations as the response (back reaction) effected by the gravity waves. 

Both the gravity waves and their back reaction on the metric field can 
be parameterized in an infinity of ways, refiecting the symmetry of the exact 
theory under generic gauge transformations. Consistent with the expansion 
above, we write gauge trasformations at first and second order in the form 

^9t^u = 6 g gf} , (2) 
9?^ = 9jS - ^iiW^9f.u + -^Iwgj^J - -^^(2)g^^^ , (3) 

where C is the Lie derivative and ^(^^'^ are two independent vectors. 

Their significance becomes clearer if we regard this gauge transformation as 
being generated by a second order coordinate transformation: 

,2 

In practical terms, the statement is that the parametrization of the first- 
order quantities (the gravitational waves) can be carried independently of 
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the parametrization of the second-order quantities (the back reaction) 

Back reaction is the feedback effect driven by the nonhnearities of the 
Einstein field equations. At second order in the perturbations, we expand 
Einstein's equations as 

G^^m+G^^[5g] + ... = Q^ (5) 
where the gravitational waves 5g obey the Einstein field equations in vacuum, 

G^ilM = ^- (6) 

The notation should be clear: G^^} is a differential operator which is linear 
in the arguments, and for this reason appears in Eq. (|^) acting on the 
first order metric perturbations, and in Eq. (^) acting on the second order 
metric perturbations. is a differential operator which is quadratic in its 
arguments. 

We will consider Eq. for the spherically symmetric (isotropic) second 
order perturbations. The Schwarzschild background metric in a; = (t, r, 9, cj)) 
coordinates is 

^7(°) = diag(-^,pr^r^in2e), (7) 
with r = 1 — R/r and R = 2GM. The metric to second order is given by 

9,. = g^^Xr) + e5g,.{x) + ^gf}{t, r) + 0{e') . (8) 

We have omitted all second-order anisotropic terms from this expression be- 
cause at this order in perturbation theory they decouple from the isotropic 
terms. 

Note that by keeping the spherical symmetry in this problem does not 
imply the time independence of the second order metric, which includes ef- 
fects from back-reaction: the time independence of the Schwarzschild metric 
is the consequence of Birkhoff 's theorem, which assumes no source terms for 
the Einstein equations. Therefore, the metric of a black hole that radiates 
can no longer be static. 

^Note, however, that the parametrization at the first order does affect the second-order 
quantities through nonlinear terms - e.g., the middle term on the right hand side of (^). 
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The back-reaction equation we are interested in is the isotropic projection 
of Eq. @, which we recast as: 

GfAg^'Kt^r)] = -(G(f)N)n = SttGT^J^ , (9) 

where the angle average is given by 

{F)n = ^Jded(j)F. (10) 

We call T^fJ^ the effective energy-momentum tensor (EEMT) of gravitational 
waves in the Schwarzschild background. 

The gauge transformation at the second order (|^) is rewritten as: 

g^^{t,r) = gf}{t,r) - {C^,.,5g,^)^ - \c^,.,gf} + \{Cl,,g''^})n . (11) 

where now ^^^^ = [^'f'\t, r), f^^)(t, r), 0, 0], i.e. only the gauge freedom accord- 
ing to the spherical symmetry remains. It will be clear in the next sections, 
when we explicitly consider the expansion of the perturbations in spherical 
harmonics, that Eqs. (|^) and (0) follow from the averaging of Eqs. (|^) 
and (^ respectively - since the second order quantities appear linearly in 
the second order relations, only the monopole part of g'^^}} and S^jf' survive to 
the averaging, and these are precisely g'^l}}{t^r) and The monopole part 
of the metric perturbation contains four degrees of freedom, but we can use 
the two components of ^^^^ to cancel two of the degrees of freedom, so we 
are left with two independent degrees of freedom. This is the right number 
of functions needed to describe a time-dependent isotropic metric, as for the 
case of an evaporating black hole 0. 

Finally, let us comment on three basic points. 

First, it is easy to verify that the EEMT obeys conservation equations 
which are just the Bianchi identities expanded to second order in perturbation 
theory. These conservation equations are completely independent of the form 
of the second-order degrees of freedom, since the part of the Bianchi identities 
which is linear in g'j}} and its derivatives is zero identically (see also 0). 

Second, we do not address the problem of the gauge dependence (or 
independence) of the EEMT of gravitational waves. Once one has fixed 
his gauge choices at the first and second order, one can then try to extract 
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physical information from the analysis of the second order Einstein equations 
or from considering gauge-invariant observables. 

Third, in what follows we consider the problem of linearized perturbations 
in vacuum. Classically, we would not have any physical solutions without 
source terms. Here we could think of the gravitational waves as having been 
generated by quantum effects close to the horizon (Hawking evaporation p]]). 

3 Linearized perturbations in vacuum 

Consider the Einstein field equations in vacuum, linearized around a back- 
ground geometry gj^J: 

+5g,M'^-g^^5g^^R'^'^^^ = Q (12) 

where 5g is the trace of dg^,^ and | denotes the covariant derivative with 
respect to the background metric. By subtracting from the above equa- 
tion its trace and considering Ricci flat background spacetimes (such as 
Schwarzschild) we have 

•^Vl"'" ~ ^%c.\u" " ^3ua\^'' + ^9\u^ = , (13) 

and 

59j'''-5g;: = Q. (14) 

It is easy to see that the last equation is trivial if one uses the transverse 
(5(7^^'" = 0) and traceless {6g = 0) gauge. By using the relation 

^Vl J" = + + Ri'^Sg,^ (15) 

we can switch the order of derivatives in Eq. (p!3|): 

'^V I J" - V - ^9.^ 1; + ^R''f!^'^S9a, + 5g\., = 0. (16) 
We can associate a scalar product with these equations of motion p!0| : 
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= -^ / d'x^n^ (17) 



s 



X 



I" 



where if) and are two generic (complex) solutions to the equations of motion 
(|T6|), is the unit vector normal to the spacelike hypersurface S, is the 
determinant of the induced metric on E and ipp^y is defined by 

i'.u = i^,u - . (18) 

The normalization of a gravity wave mode should be given by the inner 
product ([T7|). However, because metric perturbations diverge at infinity in 
the RW gauge, it is common practice to neglect the Klein-Gordon inner 
product (|T7|) and instead to normalize the gravity waves through their mass- 
energy in an asymptotically flat region of space. Since we would like to 
compute the mass-energy of the gravitational waves from first principles, 
we will avoid this procedure. Instead, we show how to obtain asymptotic 
expressions for the gravitational waves in Regge- Wheeler gauge, and use 
them in the inner product ([T7|), which then becomes well-defined. 



4 The Regge- Wheeler gauge 

Under a gauge transformation (^, the metric perturbations are transformed 
according to the general laws (||) and (H). In order to choose a gauge we must 
specify constraints for the metric components that uniquely fix the eight free 
functions ^^^^^ and ^(^^^ with respect to an arbitrary gauge. The penalty for 
not fixing the gauge completely is having to deal with unphysical degrees of 
freedom (which become ghosts after quantization), corresponding to the free 
functions that were left unconstrained. 

To first order in perturbation theory, one of the (infinitely many) choices 
of coordinates which fixes the gauge completely in the Schwarzschild back- 
ground is the Regge- Wheeler gauge [|Tl| , |], ||. This choice of gauge is excep- 
tionally convenient, since the tensor structure of the Einstein field equations 
decouple from their angular dependence. The two degrees of freedom of 
gravitons in vacuum correspond to two (orthogonal) sets of perturbations. 
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odd and even (or electric and magnetic), depending on how they transform 
under parity {x — >• —x). 

In the spherically symmetric background (0), the perturbations are ex- 
panded in spherical harmonics Y^{6, 0) and therefore the modes carry mul- 
tipole numbers {i,m). In Regge- Wheeler gauge we have, for the odd metric 
perturbations: 



^0 



sym sym 
y sym sym 

where /ig™ ^'^d hf^ are functions of t and r. 
The even parity perturbations are given by 



rig 


dYiT 


sin 9 


d<t> 
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(19) 
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Yl 





















(20) 



where i^Q™, L^"*, and K^"^ are also functions of t and r. 

With the metric perturbations written in this form, the equations of mo- 
tion can be separated and cast in the simple form (see for example |TT], ^ |^) 



Q2^o,e Ql^o^i 







(21) 



where the superscripts o, e correspond to the odd or even perturbations and 
r* is the so-called tortoise coordinate, 



r 



r* = r + R\\i\ — 
\R 

The odd and even potentials V°'^ are given respectively by: 



Vir) 



r 



3i? 



(22) 



(23) 



and 
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_ 2A^(A + l)r3 + 3A^i^r^ + |Ai^V+|i?3 

where A = £(£+l)/2-l. 

Both the even and odd potentials have the same asymptotic form in the 
hmit r — > +00: 

V°'\t)^2^^. (25) 
Both also go to zero on the horizon, but with a different slope: 
, ^2A-1 ,,,, , 2rA2 + A + 3/4 

The two odd functions /iq™ and /i^™ are easily expressed in terms of the 
Regge- Wheeler function x°'- 



_9_ 

dr* 
r 



(27) 



The even metric perturbations can be written in terms of the Zerilli func- 
tion x'^ ^ well. Since one of the linear equations of motion states that 
— = H^'^, there are only three different nonzero metric compo- 
nents: 



Kim 
Him 



dx^ 
dr* 



d_ 
di 
d_ 
dr 



9e{r)x^ + 



rge{r)x' + rr 



r dx'^ 
r dr* 
dx^ 
dr 



(28) 



where 
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_ A(A + l)r^ + 3/2i?rA + 3/2/?^ 

= r^{Xr + 3/2R) ' ^^^^ 

1 Xr^ -3/2XRr -3/4R^ 
^ ;r Xr^3/2R ■ (»' 

Some remarks are in order. Note that the potentials V°''^{r) vanish asymp- 
totically both on the horizon and at infinity (r* — =Foo). On these asymp- 
totic regions the solutions to the wave equation (0) can be expanded in 
terms of plane waves in retarded or advanced time, 

where the plus holds for ingoing (towards the horizon) waves and the minus 
holds for outgoing waves. 

This implies that in Regge- Wheeler gauge the odd and even metric pertur- 
bations diverge both on the horizon and at infinity: for example, hi ~ r*x° 
when r* — i> CX3 and hi ~ Re~'^*^^x" when r* — oo. Of course, this is a co- 
ordinate artifact of the Regge- Wheeler gauge, and physical quantities should 
remain finite in that (or any other) gauge, as long as we keep away from 
the singularity at the center of the black hole. There are gauges in which 
physical quantities are manifestly regular at infinity, such as asymptotically 
fiat P] or radiation gauge f^. Our purpose is to show (see Section 5) that 
the Regge- Wheeler gauge can also be used to compute physical quantities, 
even if the metric components diverge. 

We can simplify greatly the ansatz for the metric if we adopt the standard 
procedure of rotating the 2;-axis to put each mode (£, m) in the state (£, 0). 
The dependence thus drops out of the metric perturbations (|TW)-(PD|), and 



all subsequent equations involve only Legendre polynomials Pf(cos6') and 
their derivatives. Of course this simplification does not change our results 
for the back reaction on the isotropic mode. 

If the metric ansatz (|l^)-(|3) is appropriate to describe the anisotropic 
metric perturbations around the black hole, what is appropriate to describe 
the back reaction of these perturbations on the isotropic background? Physics 
dictates that we should expect some mass to be lost by the black hole. 
However, the mass appears both in qoq = —(1 — 2GM/r) and in gu = 
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(1 — 2GM /r)^^ . Generically there should be two degrees of freedom describ- 
ing the (no longer static) isotropic metric||], which we define adopting the 
following gauge choice at the second order: 

-1 

+ 0(e3) . (33) 

The interpretation of AqM is related to the ADM mass at infinity, while 
AiM is related to the location of the horizon^. However, sometimes it is 
useful to make the following approximation: assume that the black hole 
emits radiation in packets, such that for a spherical shell of radius r = f 
from the black hole, any given packet is either completely inside or outside 
the shell. If the gravity waves packet is inside the shell, the metric describing 
the isotropic background is just the exterior Schwarzschild metric without 
corrections. If the packet is outside the shell, the mass of the black hole 
measured by an observer located at r = f has changed by an amount AM 
that is now independent of time and radius. In the limit of this approximation 
we can ignore subtleties with the dynamics that may distinguish AoM(r, t) 
from AiM(r, t) in a more complex physical situation. 

Let us comment on a peculiarity of the gauge choice at second order. 
There is no way of fixing completely the gauge for the isotropic perturbation, 
since the remaining freedom is related to a redefinition of time: 

t^t + p{t) + 0{e') (34) 

This is also related to the fact that time is defined ad hoc for the Schwarzschild 
solution (which for us is the zero order metric of our perturbative analysis). 



9ii 



2GM ^2GAiM{r,t) 



+ e 



5 The flux of gravitational waves 

In this section we will focus on the mixed t — r component of the EEMT of 
gravitational waves for two main reasons. The first is simplicity: G*^ = R^^ 

^We thank Roberto Casadio for having pointed out this feature to us. 
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to second order has a short expression in terms of the perturbations. The 
second is that this term represents an energy flux. Consider for example a 
scalar field in the Schwarzschild background. The mixed t — r component of 
the scalar's energy-momentum tensor is given by: 



T* = ^<^,t <^,r* (2^ - 1) + 2^(/? (^ip,tr* -^y^,t ) , (35) 

where ^ is the coupling of the scalar field to the curvature. By considering 
the proper normalization ~ 1/r in front of the asymptotic spherical plane 
waves solution for a scalar field one has an asymptotic behaviour (9(l/r^) on 
the horizon and (9(l/r^) at infinity. We observe that the behaviour of on 
the horizon is not really pathological, indeed it is regular in a freely falling 
frame. 

We now proceed to calculating the mixed t — r component of the effective 
energy-momentum tensor (EEMT) of gravitational perturbations around a 
black hole. Since the even and the odd degrees of freedom are orthogonal 
to each other at this order in perturbation theory, we can consider their 
contributions to the EEMT separately. For the choice of second order metric 
coefficient in Eq. ( p3D one has: 

R'}'V\t,r)] = -'-^ (36) 

The odd contribution to the t — r component in Regge- Wheeler gauge is, 
after averaging over angles: 



E 



i{i + 1) /ii 

2^+1 2r4 



r'^{h\ - hi) - 2rhi + 2/io 



(37) 



where a dot and a prime denote derivatives with respect to t and r respec- 
tively. 

The first term inside square brackets in (|37|) appears to go as at r — oo 
(since ho, hi oc r in this limit). If that was so, the radiated mass AM would 
be divergent. In reality, that term is at most oc because the leading term 
cancels upon use of relations ([19| ) and the equations of motion (^). Nev- 
ertheless, the C(l/r) term in {R\.odd)^ still gives a contribution that could 
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make AM divergent. We show in the next section that these apparently 
divergent terms in Regge- Wheeler gauge are similar to the terms in the nor- 
malization of the gravity waves that are also apparently divergent. However, 
in both cases the divergences are revealed to be fictions, and the physical 
information can be retrieved. 

On the horizon the leading term in the right hand side in ( pTf ) appears to 
go as 0{l/r^). Also in this case by using the equation of motion (|2T|) this 
term vanishes, and so T^/^^ ~ (9(l/r^), as in the scalar field case (pS)). 

The even contribution to the t — r mixed component has a similar struc- 
ture: 

{Rle.en)n = ^^^^ [^^' ^' ' H) + 2r'T K" L (38) 
+r^K'{k + H)- r^KH' + ArVK'L 

Again, from ( ^8]) there is an apparent divergence in T\!^^ at infinity which is 
given by terms 0{r^) inside the square brackets in (|5BD. However, the leading 
term vanishes by using (^Tf ), and the physical information can be obtained 
among the subleading terms. Also in this case the apparent divergence on 
the horizon 0{l/r^) vanishes by using (^T]) . ^ 



6 Asymptotics of the Regge- Wheeler and Zer- 
iUi functions 

Consider now the normalization of the metric perturbations in Regge- Wheeler 
gauge. Since perturbations with different parity decouple the odd contribu- 
tion to the inner product (p!7| ) can be written as: 



{r',MoM = -^E^^^//^*r (39) 

^ Note that all these cancellations of the most divergent terms by using the equations 
of motions are also a property of the energy momentum tensor of gravitational waves, even 
before the average over angles. 
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X 



hihl - hlhi + ^{rh[ + 2hi) - y (r/if + 2hl] 



where the angle integrations have aheady been performed. When r — > oo 
the integrand appears to be oc r^. This structure appears in the the even 
contribution to the inner product as well. 

We show next that in fact no divergences survive in the normalization 
conditions once the perturbative solutions to the wave equation (|21| ) are 
used in the expression above. 

From Eq. (^) for the Zerilli {x^) and Regge- Wheeler functions it is 
clear that at infinity they assume the form of a plane wave in retarded time 
u = t — r* (or advanced time f, in the case of an incoming mode). This 
suggest an expansion of x terms of functions of both retarded time and 
radius of the form^] 

X(r, t) = Xo{u) + r-'X^{u) + r-^X^iu) + 0{r-^) . (40) 

This expansion should be consistent as long as the gravity wave packet does 
not probe the region where the potential is large, r ~ 3GM. In other words, 
the approximation breaks down as soon as the reflection and transmission 
of the waves become important and the very description of a wave packet in 
terms of purely ingoing (or outgoing) modes breaks down. 

We can substitute the ansatz above into the wave equation (|2l|) and solve 
the hierarchy of equations that ensue. Since the leading order terms in the 
potentials V° and are the same to order r~^, the solutions for x'^ and x" 
are identical to each other at that order. The result is, for an ingoing gravity 
wave of even or odd parity. 



X°'1r,t) = ^X^'0^u) + r-'X^'''{u) (41) 



2 ^ ^ 2(A + 1) ^ ' 



+ 0{r- 



Notice that X is a completely generic function of retarded time, not nec- 
essarily a plane wave, reflecting the fact that we have not yet imposed any 
boundary conditions on the problem. We also stress that X carries an index 
i, like the mode of metric perturbations. 
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From now on we focus on the odd perturbations for the sake of simphcity, 
since the expressions for the even perturbations are too cumbersome, and the 
physics and the lessons we draw are precisely the same as in the odd case. 



With the help of relations ( pTf ) we can express the metric perturbations in 
terms of the odd function X{u): 



ho 

hi 



A + 1 

r 



-X 



A 



A + 1 



-X 



-1 



GM 



2(A + i; 



X - Ax 1 + 0(r-2) ,(42) 



A + 1 

-1 



+r 



A + 1 2(A+1) 2 ^ 



(43) 



C(r 



-2\ 



It is now a matter of algebra to substitute these expressions into Eqs. 
( PTD and (psp. The result for the normalization is 



KA 



2i+l 



dr* 



A d 
'{X + Tydi 



(x^x;) + 0(r°) 



(44) 



Therefore, the 0{r^) divergence in the integral (|39D has simply disappeared 
after we used the asymptotic formula (p]). 

The first term in the expression above would make the integral diverge 
as r^. However, that term vanishes because it is the time derivative of an 
integral over the spacelike hypersurface S (our assumptions are that the grav- 
itational wave packet is either entirely inside or entirely outside S.) There is 
a physically intuitive way to argue that this divergent term should vanish. In 
the WKB approximation the X^'s are plane waves, Xa oc exp{—iwau). In this 
framework, the phases in the integrand of (|4^ ) would interfere destructively. 

In either case, the normalization of the odd metric perturbations is given 
by the next term in perturbation theory: 



24A2 



a/3/odd,WKB 



(A + l)(2£+l) 



dr*ix^x: 



(45) 



If the mode solutions X^ and X^ are just WKB modes, the integral above 
gives the expected normalization ||10| with a delta function over the frequen- 
cies of the modes, 6{w^ — w^). 
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Therefore, it is possible to normalize gravitational waves in Scliwarzscliild 
background using the inner product in Regge- Wheeler gauge. As we show 
below, a similarly divergent time derivative term appears in the equation for 
the time variation of the mass of the black hole. Again, the next to leading 
term is the finite one, in agreement with a finite gravitational radiation flux 
from an object with finite mass. We note that the term in ( ^41) which is 
proportional to r at infinity contains also a divergence on the horizon: we 
have just shown how deal with it. The next term in (^5]) is regular also on 
the horizon, since an expansion in powers of F similar to (^) holds close to 
the horizon. 

Let us consider now the flux of energy given by Eq. (|37|) . Substituting 
the expressions (^) and ( |43| ) we obtain the following (apparently) divergent 
term: 

A.M^I: (,^1)P,^1) 4a-A- + 0(.°). (46) 

Some authors have assumed [P, Q that X obeys some asymptotic condi- 
tions such that the 0{r) term vanishes. It appears to us that the WKB 
approximation could already be sufficient to ensure the finiteness of the ra- 
diated mass. 

The radiated mass from odd-parity gravitational waves comes from the 
next order term, and using (|4^)-(|43|) the final result is 

A.M-i:^/^*(rf + 0(.-), (47) 



which is the expression that has been known for a long time |T2| and was 
first found using the Landau-Lifshitz pseudo-tensor in the asymptotically 
fiat gauge. Other authors 0, Q have obtained similar results for the even 
perturbations (in the 1 = 2 multipole case), also using the Regge- Wheeler 
gauge. 



7 Conclusions 

We have explored the possibility of performing back-reaction computations 
for the Schwarzschild black-hole in the Regge- Wheeler gauge. The Regge- 
Wheeler choice has the advantage of fixing completely the gauge, but it has 
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the disadvantage of divergencies in the metric coefficients both on the horizon 
and at spatial infinity. 

We have shown that divergencies in the gravitational energy fiux at infin- 
ity and on the horizon, which could arise from choosing the Regge- Wheeler 
gauge, are in fact nonexistent. Furthermore, the same types of cancellations 
arise in the evaluation of the inner product, which is related to the equations 
of motion of the linear perturbations. This could represent a viable check 
for the amplitude of gravitational waves, independent from the requirement 
that the energy fiux of gravitational waves must be finite at infinity. The 
investigation of the effect of the back-reaction of gravitational radiation on 
the horizon is under current investigation. 
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